The Spin-Exchange Dynamical Structure Factor of the 5 = 1/2 Heisenberg Chain 
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We determine the spin-exchange dynamical structure factor of the Heisenberg spin chain, as is 
measured by indirect Resonant Inelastic X-ray Scattering (RIXS). We find that two-spin RIXS 
excitations nearly entirely fractionalize into two spmon states. These share the same continuum 
lower bound as single-spin neutron scattering excitations, even if the relevant final states belong 
to orthogonal symmetry sectors. The RIXS spectral weight is mainly carried by higher-energy 
excitations, and is beyond the reach of the low-energy effective theories of Luttinger liquid type. 



A remarkable feature of quantum systems in one di- 
mension (ID) is that quantum criticality is the norm 
rather than the exception pQ. ID and quasi-lD systems 
as diverse as carbon nanotubes, stripes in cuprate high- 
temperature superconductors, confined ultracold atomic 
gases and quantum spin chains all provide realizations of 
critical quantum liquids. Possibly the most studied pro- 
totypical ID quantum critical system is the Heisenberg 
5 = 1/2 antiferromagnetic chain [3J. Its basic excitations 
are spinons, fractionalized spin excitations that emerge 
in the critical state [3]. Until recently, the only tech- 
nique available to investigate these fractional spin exci- 
tations was inelastic neutron scattering (INS). In neutron 
scattering integer spin-flip excitations carrying S—l are 
created, so that the INS response involves the pairwise 
creation of spinons. The INS amplitude is determined 
by the single-spin dynamical structure factor (DSF) and 
a precise matching between the observed INS intensity 
and the dynamical structure factor calculated from the- 
ory has recently been achieved gHB] . 

A few years ago, x-ray photon scattering emerged as a 
new tool to measure dispersive magnetic excitations [TJ. 
When the incident photon energy is resonant with an ab- 
sorption edge of the material, yielding Resonant Inelastic 
X-ray Scattering (RIXS), the correlations between two 
neighboring spins factor into the magnetic x-ray scatter- 
ing processes. The corresponding scattering amplitude 
is given by the momentum dependent two-spin dynami- 
cal structure factor, also referred to as the spin-exchange 
DSF. When for instance in an antiferromagnetic copper- 
oxide compound the incident energy is tuned to the cop- 
per K-edge - so-called indirect RIXS - exclusively this 
type of magnetic scattering occurs [THTT]. The recent 
experimental advances in inelastic x-ray scattering allow 
to probe one of the most elementary, fractionalized states 
of critical matter in an entirely new manner, highlighting 
more elaborate spin-spin correlations. 

It is the purpose of this Letter to provide a nonpertur- 
bative calculation of the spin-exchange RIXS scattering 
amplitude for the Heisenberg 5 = 1/2 antiferromagnetic 
chain 

H = J^(S,-S i+1 -l/4). (1) 



Our results aim on the one hand to challenge, guide 
and inspire further experimental magnetic RIXS efforts; 
on the other hand the response function that we con- 
sider here is of fundamental significance and, as explic- 
itly shown later on, beyond the reach of traditional low- 
energy based theories (such as e.g. Luttinger liquid the- 
ory) because it sits at energies of order of the exchange 
J. Our problem and approach thus represent an uncom- 
mon example of direct contact between experiment and 
theory beyond universality. 

The model's integrability [T2TH5] is well-known to per- 
mit nonperturbative calculations of equilibrium proper- 
ties, but it has recently also become possible to accurately 
compute certain dynamical properties for both finite |16l - 
|2"0] and infinite [2TH21] systems. Until now these predic- 
tions have been limited to the single-spin DSF measured 
by INS. For experiments other than INS, governed by 
different correlators, at present no results are available. 
Here we fill this gap for RIXS, using the exact Algebraic 
Bethe Ansatz to calculate the spin-exchange DSF. 

Magnetic RIXS cross section - Before calculating its 
response function, we briefly describe the magnetic RIXS 
process. RIXS is a photon-in photon-out scattering tech- 
nique in which the energy of the incident x-ray photons is 
tuned to an atomic absorption edge of the material that is 
studied. X-rays can for instance be tuned to the -FT-edge 
of a transition metal ion (1 — lOkeV), for example copper, 
in the sample. In this case a magnetic scattering process 
as sketched in Fig. [Tjcan occur [E1[H]. An x-ray incoming 
on site j produces a Is — Ap electronic transition and cre- 
ates a core-hole on the jth copper ion of the chain. In the 
intermediate state, the presence of the core-hole modifies 
the 3d on-site energy levels through the Coulomb inter- 
action. In the Mott-insulating limit, this perturbation 
modifies the spin-exchange process with the two neigh- 
boring 3d electrons, and thus locally modifies the actual 
superexchange coupling J between the spins. Following 
the notation of [H E] we denote the perturbed coupling 
J c = (1 + rj) J. In the final state, the core- hole is filled 
again by the Ap electron, but the spin chain is left behind 
in an excited state. The cross section for this scattering 
process is given by the Kramers-Heisenberg relation. As 
a function of energy loss to = — ujj ut and momentum 
transfer along the spin chain q — qi n — q out the scat- 
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FIG. 1: Mechanism by which double-spin flip transitions are 
created in the indirect magnetic RIXS process. 
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Here \n), \f) are the initial, intermediate and fi- 
nal states with respective energies Ei, E n , Ef and D 
is the dipole operator that creates or annihilates photo- 
excitations. Because the Is core-hole is highly energetic, 
it quickly decays, leading to an energy broadening T of 
the intermediate state - its inverse lifetime. For T ^> E n , 
A fi can be expanded in a power series which re-summed 
to leading-order ultimately provides the x-ray scattering 
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S exc (q,u>). It incorporates 
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the spin-exchange dynamical structure factor 



S exch (q,u) = 2tt' 



X q \a)\ 2 5{u-Lo a ), (2) 



where the ground state is |0), the excited states \a), ex- 
citation energies uj a — E a — Eq$ and the spin-exchange 
operator is X q = £\ e i<?i ( S i-i ' S j + S j ' s i+i)- Mo- 
tivated by the successful correspondence between pre- 
dictions and experiments for the long-range ordered 2D 
Heisenberg antiferromagnet [9-1TH2S], we se ^ out to com- 
pute this structure factor in the quantum critical ID case, 
applicable for instance to Sr2Cu03. 

Computing the spin- exchange DSF - So far exact cal- 
culations for the spin chain have been restricted to the 
single-spin DSF, as measured in neutron scattering: 



|S»| 2 5(a; 



■Wa)i (3) 



where (0|S^|a) is the form factor (FF) of the Fourier 

transformed spin operators S q — (1/V^V) £,-=i er^S®, 
a — z. +, — between the ground state and eigenstates of 
excitation energy uj a (see Fig. [2]). The strategy which 
we adopt for the calculation of the RIXS response func- 
tion S exch (Eq. [| is that of the ABACUS method 
[26] • Within this approach eigenstates are explicitly ob- 
tained from Bethe Ansatz, matrix elements from the Al- 
gebraic Bethe Ansatz, and the trace over intermediate 
states is performed numerically using an optimized search 
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FIG. 2: Single spin dynamical structure factor S"T gle (q,uj) 
of the Heisenberg chain. The DSF is dominated by the 2- 
spinon spectrum but also 4-spinon and higher excited states 
contribute (6%). The computation is for N = 400 sites. 



through the Hilbert space. In what follows, we briefly de- 
scribe each of these three steps for the specific case of the 
RIXS response function (Eq. [2]) . 

The Hilbert space can be divided into subspaces of 
fixed magnetization characterized by the number of 
downturned spins M. Eigenstates of Eq. [I] for an (even) 
N sites periodic spin chain are completely determined for 
M < -j by a set of M rapidities {Ai, . . . , Am} which for 
the isotropic chain solve the Bethe equations 

7T 1 M 

arctan (2A,) = —U + — ^ ar ctan (Ai - A fe ) (4) 



k=l 



with I\ , . . . , Im a set of integers for odd M and half- 
odd integers for even M. Each set of quantum numbers 
specifies a set of rapidities. The ground state is defined 
by {J° = fc - ^±1 }, k = 1, . . . , M. The energy and the 
momentum of an eigenstate are given by 
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Turning now to the matrix elements of the exchange 
operators X q , we can exploit the spin isotropy of the 
system to express the full spin-exchange DSF matrix el- 
ement in Eq. [2] as a function of the SfS* +l FFs only. 
By globally rotating the (0| SfSf +1 \a) and (0| S]S] +1 \a) 
FFs appearing in Eq. [2] about the y and x spin axes, re- 
spectively, and using the fact that the ground state is 
a global su(2) singlet, one can show that only singlet 
excited states contribute to Eq. [2j This conclusion is 
reached by first noticing that the SfSf +l operator creates 
excited state only in the S to t = 0, 1, 2 sectors. Secondly, 
the rotation of states belonging to these sectors and with 
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FIG. 3: Spin-exchange dynamical structure factor S exch of 
the isotropic Heisenberg chain. The sharp on-set of the spinon 
continuum is very similar to the single-spin DSF S^ 9 . The 
intensity around q — tv is suppressed, thereby enhancing high- 
energy response around q — f , Tp. The computation is for 
N = 400 sites. 

zero magnetization (S% ot = 0) gives 

s " \a, Sf ot = 0, Sua = 0) = \a, 0, 0) 
e^i^ \a,0,l) = cos(0)|a,O,l) + .. . 

j^v > l + 3cos(20), 
e^ s >|a,0,2) = _AZZ| a7 o,2) + ... 

(in which . . . represent Sf ot ^ states of which the 
SjSj +1 matrix elements with the ground state vanish) 
with similar results for a rotation around the y axis. With 
rotations of respectively ±7r/2 to transform either 5| or 
S? into Sj, the total contribution of eigenstates with 
Stot = 1)2 vanishes, allowing one to rewrite the spin- 
exchange DSF as 

^(q^)=cos 2 (g/2)^x 

£ ^|e^"(0|^ +1 |a)| 2 <5( W - WQ ). (5) 

aeS tot =0 3 

In order to calculate the matrix elements (0| S*S* +1 \a), 
we make use of Algebraic Bethe Ansatz methods [TJl |2~T1 - 
1317] to represent these as explicit functions (taking the 
form of matrix determinants) of the rapidities involved 
in the left and right eigenstates. State norms are also 
given by an explicit determinant |13|. 1311 152") . For brevity, 
explicit expressions are not presented here |33j . 

The third step, summation over intermediate states, 
can now be performed to obtain quantitative results 
for the DSF in Eq. [2] This is done using the ABA- 
CUS algorithm [26] which sums intermediate state 
contributions in a close to optimal order. By solving 
the Bethe equations, the eigenstate rapidities give 
the energy, momentum and the FF value through the 



FIG. 4: Fixed momentum profiles of the spin-exchange dy- 
namical structure factor S exch (q,uj)/(cos 2 (q/2)) (plain) and 
S°l ngle [q ,0;) (dashed), each normalized to its own sum rule. 
Data are plotted as function of the energy loss oj, at q = 7r/4, 
tt/2, 3tt/4. 



determinant expression. Large families of excited states 
are summed over until satisfactory saturations of sum 
rules are achieved. More precisely, first, we compare the 
computed integrated intensity with the analytical result 
in the infinite chain limit: / & i Y, q S exch {q, u) = 
| — ln(2) + |C(3). A second check is the first 
frequency moment sum rule J ^u}S exch (q,uj) — 
6sin 2 (g) {(a* - x 2 ) (l - 4cos 2 (g/2)) + 3C(3)-4in(2) j 

(valid for the isotropic case, in which Xi = (SjSj+i); 
we have made use of the results of [34]). To obtain 
smooth curves in frequency ui, the delta function in 
Eq. [2] is broadened to a scale commensurate with the 
energy level spacing. The results for the RIXS scattering 
amplitude are presented in Fig. [3] (for all momenta in 
the full Brillouin zone) and [4] (as a function of energy at 
fixed momentum). 

As the RIXS response involves two spin excitations, to 
each of which a S = 1 is associated, one might naively 
expect it to be dominated by four spinon intermediate 
states. However, it turns out that the RIXS excitation 
fractionalizes almost completely into two spinons. The 
state contributions to the sum-rule shows that the two- 
spinon states which cover all but 10 _4 % of the RIXS 
signal for a finite size chain with N = 400. This ra- 
tio is even more drastic than for the INS results where 
the excitation of four and more spinons is responsible for 
6% (N = 400) of the signal. Interestingly, the magnetic 
RIXS response is at first glance rather similar to the neu- 
tron scattering one, sharing the same lower bound of the 
spinon continuum, even if the RIXS and INS final states 
belong to different symmetry sectors and are orthogo- 
nal. The fact that the spin-exchange DSF exists within 
a continuum is a direct and explicit consequence of the 
fractionalization of spin excitations in the quantum crit- 
ical spin chain. Even if the excitation continuum probed 
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by RIXS exactly coincides with the one probed by INS, 
the spectral weight has a markedly different distribution. 
This crucial difference is partly caused by the static fac- 
tor cos 2 (g/2) which originates from the modification of 
two neighboring exchange couplings in the x-ray scatter- 
ing process. Excitations are thus associated to a typical 
length 2a (with a the lattice spacing) or equivalently pre- 
dominantly carry ±5 (mod 2tt) momentum. The figures 
clearly illustrate this fact: for INS (Fig. [2]) the signal is 
maximal at the antiferromagnetic wavevector q — ir at 
lo close to zero, whereas the RIXS amplitude vanishes 
there. Rather it is concentrated above the continuum 
threshold at q = tt/2, characterized by a vanishing group 
velocity. The fixed-momentum profiles in Fig. [4] show 
further differences between RIXS and INS apart from 
the cos 2 (g/2) factor. While the weight of both responses 
predominantly sits between the top and bottom of the 
two-spinon continuum, the weight distribution between 
these two are clearly distinct, the RIXS response having 
a much broader shoulder at higher ui than the INS re- 
sponse. This quantitative difference is sufficiently large 
to be observable experimentally. 

The asymptotic behavior of correlations along the spin 
chain are well described by low energy effective theories 
PQ, which provide a detailed understanding of the anti- 
ferromagnetic singularity at momentum 7r and zero en- 



ergy in the INS response. The exponent of the singular- 
ity at the lower bound of the spectrum is given by the 
asymptotic decay of the two-spin correlation. Straight- 
forward application of the known results for the four-spin 
correlator from Luttinger liquid theory [35] to the RIXS 
intensity fails as the signal vanishes at all low energies. 
The development nonlinear extensions to Luttinger liquid 
theory [36 or alternately the quantum group approach 
[2"TH2"4"| might bring this within reach in the future. 

We have, in conclusion, explicitly calculated the rel- 
evant magnetic response function for indirect RIXS on 
the Heisenberg chain, which involves the SiSj+i spin- 
exchange operators. This operator predominantly probes 
high-energy magnetic excitations and the integrability- 
based method which we have employed here is then 
the only one capable of describing RIXS. Alternate ap- 
proaches based on approximate low-energy theories can- 
not reach energies of order of the exchange J. We have 
demonstrated in more general terms that via the Alge- 
braic Bethe Ansatz formalism a new family of correlators 
of the type SfS^ +1 with a,b = z,+,— is accessible for 
computation. The challenges of future work are to ex- 
tend the present exact calculation of the spin-exchange 
dynamical structure factor to anisotropic spin chains, in- 
cluding also finite magnetic fields. 
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